Systematic Study of Gravitational Waves from Galaxy Merger 
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A systematic study of gravitational waves from galaxy mergers, through iV-body simulations, 
was performed. In particular, we investigated the relative importance of galaxy components (disk, 
bulge and halo) and effects of initial relative velocity, relative angular momentum and mass ratio 
of the galaxies. We found that the features of light curve of gravitational waves, such as peak 
width and luminosity, are reliably simulated with particle numbers larger than ~ 10 4 . Dominant 
contribution to gravitational wave emission came from the halo component, while peak luminosity 
amounted to 10 31 erg/sec for the collision of two halos with mass 3.8 x 10 la /i _1 M o . We also found 
that the initial relative velocity in the direction of the initial separation did not significantly affect 
gravitational wave emission, while the initial relative angular momentum broadened the peak width 
and suppressed the luminosity. Mass dependence of the peak luminosity was also investigated, and 
we obtained evidence that the luminosity is proportional to the cubic mass when the scaling relation 
is satisfied. This behavior was considered by a simple analysis. 

PACS numbers: 



I. INTRODUCTION 

Structure formation of the universe proceeds through 
gravitational interaction from tiny density fluctuations. 
In this process, gravitational waves with cosmological 
scales are expected to be produced. In [H-fl], gener- 
ation of gravitational waves was studied as a second- 
order effect of cosmological perturbations. It was shown 
that the density of gravitational wave background, pew, 
amounts to pcw/Pc ~ 10~ 20 — 10~ 15 for a wide fre- 
quency range, where p c is the critical density. At smaller 
scales, gravitational waves, radiated through the forma- 
tion of dark matter halo, were studied in Q and the 
density was estimated as pgw/Pc ~ 1CU 20 at frequen- 
cies 1CU 18 — 1CU 17 Hz. Thus, these processes of structure 
formation are imprinted on the gravitational wave back- 
ground. 

Although such gravitational wave background of cos- 
mological scales is inaccessible to direct detection, it 
could be probed by measuring the B-mode polariza- 
tion of the cosmic microwave background. In actu- 
ality, several missions suited for this aim have been 
plann ed; for example, L itcBIRD, QUIET, POLARBeaR 
(sec, http://cmbpol.kek.jp/indcx-c.html), and Cosmic 
Inflation Probe (see, http: / /www.cfa.harvard.edu/cip/| . 



Their primary purpose is to detect primordial gravita- 
tional waves generated during inflation. However, it was 
reported in [J] that gravitational waves from halo forma- 
tion could be dominant if the energy scale of inflation is 
below - 10 15 GeV. 

In this paper, we investigate galaxy merger as another 
process which produces gravitational waves of cosmolog- 
ical scales. In the hierarchical model of structure for- 
mation, it is argued that low-mass dark matter halos 
repeatedly merge with each other to form more mas- 



sive halos (see, e.g., p| @). Because halo merger is 
a highly nonlinear process which involves huge masses, 
it is expected that substantial amounts of gravitational 
waves are emitted. In Q, gravitational waves from 
galaxy merger are calculated with A^-body simulations. 
These researchers studied several typical configurations 
of galaxy merger and estimated the luminosity to be of 
order 10 33 erg sec -1 . 

By extending this previous study Q , we perform a sys- 
tematic study of gravitational waves from galaxy mergers 
through A"-body simulations with Gadget2 d, Q . While 
in 0] initial condition was generated by an open code 
"GalactlCS" and Jaffe model [ll|, here we use a 
halo adopted from a cosmological simulation as well as 
ones generated by "GalactlCS" . The former has a realis- 
tic density structure which is well approximated by the 
Navarro- Frenk- White (NFW) profile |H- On the other 
hand, the latter describe disk galaxies with three com- 
ponents, disk, bulge and halo. With these initial con- 
ditions, we investigate the relative importance of galaxy 
components (disk, bulge and halo) and the effects of ini- 
tial velocity, relative angular momentum and mass ratio 
of the halos. Also we check the dependence of the results 
on simulation resolution by varying the particle number. 
This study can be regarded as the first step toward eval- 
uating the gravitational wave background from galaxy 
merger. 

The paper is organized as follows. In section [Hi we 
explain our method of calculation, including the prepa- 
ration of the initial conditions. In section Hill the results 
are presented. Finally section IIVI is devoted to the sum- 
mary of our study and discussion including our future 
plans. 
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II. METHOD 

We follow the merging process of two galaxies using 
the parallel A-body solver Gadget2 @, II in its Tree-PM 
mode. The gravitational softening parameter is set to 
around 1/60 of the tidal radius of galaxies. We use two 
types of initial conditions (type A and B). Type A is 
prepared by generating galaxies with "GalactlCS" while 
a realistic halo is adopted from a cosmological simulation 
for type B. Detailed description of the initial condition 
will be given in the next section. 

Based on the A^-body simulation, we compute the 
amplitude and luminosity of the gravitational waves. 
Here, we neglect the effect of cosmological expansion be- 
cause the initial separation of galaxies is taken to be 
~ 1 /i _1 Mpc and the ratio of acceleration due to cosmic 
expansion and gravitational force between galaxies is of 
order 0.1. Further, the velocity dispersion of stars in each 
galaxy and the relative velocity between two galaxies are 
nonrelativistic (v/c < 10 -3 ) so that the slow-motion ap- 
proximation is valid. By denoting the deviation of the 
metric from Minkowski spacctime as and taking the 
transverse-traceless (TT) gauge, the spatial components 
of hJ^J are given by the quadrupole formula as 

^ T = f^(*-r/c), (1) 

where G is the gravitational constant and r is the distance 
from the observer to the merging galaxies. Here ly is the 
components of the traceless internal tensor given by 

j« = E"» w (W ) -|««» wa ) (2) 

(k) 

where x\ represents the position of fc-th particle. Later 
we will set the initial separation of two galaxies in the 
direction of x. 

The luminosity of gravitational waves is given by 

L<3W = (if) = I? ^ Iij [CTg S ~ 1] - (3) 

It should be noted here that this formula needs third- 
order time derivatives which may cause a substantial nu- 
merical error. We will estimate the numerical error by 
three different calculational procedures in section UlI Al 

III. RESULTS 

In this section, we present the results of our numeri- 
cal simulations. In subsection 1111 Al we generate initial 
conditions using "GalactlCS" and investigate the rela- 
tive importance of the galaxy components (halo, disk and 
bulge). Also we check the convergency of the results by 
varying the particle number. In subsection 1111 Bl we use 
initial conditions adopting a dark halo from a cosmolog- 
ical simulation. Finally we study the effects of the initial 



TABLE I: Models used in this subsection. Column(l) model 
name. Column(2):total mass. Column(3):mass of halo. 
Column(4):mass of disk. Column(5):mass of bulge. Col- 
umn(6):disk scale radius. Column(7):tidal radius. 



model 


M tot (10 1U M @ ) 


M h 


Md 


A/ b 


r d (kpc) 


n (kpc) 


model-0 


31.1 


26.5 


3.9 


0.7 


4.5 


56.5 


model-A 


31.9 


25.0 


4.4 


2.2 


4.5 


98.9 


model-D 


195.6 


188.9 


4.5 


2.2 


4.5 


330.7 



relative velocity, relative angular momentum and mass 
(ratio). 

A. Simulations with Type A Initial Conditions 

First, let us describe the initial conditions. In this 
subsection IIII Al we use initial conditions generated by 
"GalactlCS" [10| . This is based on a semi-analytic model 
for the phase-space distribution function of an axisym- 
metric disk galaxy, and outputs the distribution func- 
tions of a disc, bulge and halo from input parameters 
such as total mass and angular momentum. One of the 
advantages of using "GalactlCS" is that it can create a 
disk galaxy with halo, disk and bulge of arbitrary mass 
self-consistcntly. We will see that the dominant contribu- 
tion to gravitational wave emission comes from the halo 
component. Another advantage is that we can choose the 
number of particles arbitrarily so that we can study the 
effect of resolution on gravitational wave emission. 

We prepared three initial conditions whose parameters 
are shown in TABLE Q] Here and r t are the scale 
radius of the disk and the tidal radius, which is the out- 
ermost radius of the halo [l(| , respectively. The model-0 
is similar to the disk galaxy adopted in [7[. As we can 
see in Fig. [TJ the rotation curve is decreasing outward 
within 30 kpc, which means that the halo component 
is relatively concentrated. In this sense, this model is 
not so realistic. To study the effect of concentration on 
gravitational- wave emission, we consider model-A, which 
has almost the same mass as model-0 but a flatter ro- 
tation curve (Fig. [T]). Further, modcl-D is considered 
to show the effect of mass and halo extension. Model-D 
has a similar mass and structure as our galaxy and the 
central density profile is similar to model-A. It should be 
noted that model-A and model-D are identical to model 
A and model D in [To| . In this subsection, we will mainly 
use model-0 and, at the end, make a comparison between 
model-0 and modcl-A. 

The particle number per galaxy is set to 5.4 x 10 5 for 
all models in this section. Then we simulate the head-on 
merger of two identical galaxies with initial relative ve- 
locity in the a;-direction which corresponds to the kinetic 
energy equal to potential energy. 

Fig. [2] shows the time evolution of the luminosity of 
gravitational wave for modcl-0. We can see two major 
peaks which correspond to the collisions of two galax- 
ies. These gradually relax into a single galaxy and 
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FIG. 1: Rotation curve , V c {r) = ^GM(r)/r, for 
halo+disk+bulge (solid), halo (dashed), disk (dotted), and 
bulge (dash-dotted) for a range from galactic center to r t . 
Top: model-0. Middle: model-A. Bottom: model-D. 



gravitational-wave emission becomes much less effective 
after t = 1.2 Gyr. The peak width, ^0.1 Gyr, reflects the 
dynamical timescale of the galaxies. The peak luminosity 
reaches 1.4 x 10 30 erg/sec and the total emitted energy is 
about ~ 10 46 erg. As we will see below, the steady emis- 
sion after t = 1.2 Gyr is substantially affected by the 
simulation resolution and would be considered as noise 
within the simulation. Therefore, there is a possibility 
that more peaks do exist below this noise level. Nev- 
ertheless, it is robust that the luminosities of the third 
and later peaks are less than ~ 10 26 erg/sec and make 
negligible contributions to the total luminosity. 

As we saw in Eq. (j3]), we need the third-order deriva- 
tives of Iij to calculate gravitational wave luminosity. 
However, if we calculate these by differentiating 1^ three 
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FIG. 2: Time evolution of luminosity of gravitational waves 
for model-0. 




FIG. 3: Top: comparison of gravitational wave luminosi- 
ties calculated with three different approaches for model-0. 
See text for detail. Bottom: ratios of luminosity, F1/F3 and 
F2/F3. 



times numerically, it may cause a substantial numerical 
error. Because Gadget2 outputs well as Xi and Xi, 

it is possible to obtain I y with just one numerical differ- 
entiation, allowing for the reduction of numerical error. 
To aide in estimating the error concerning the numeri- 
cal differentiation, we show in Fig. [3] the comparison of 
gravitational wave luminosities calculated with three dif- 
ferent approaches. Fl is calculated from Iij obtained by 
Xi, F2 is calculated from 1^ obtained by Xj and x%, and, 
finally F3 is calculated from 1^ obtained by Xj, ii and Xi. 
We can see that the differences are smaller than factor 2 
during almost the entire process. Thus, we can conclud 
that the error concerning the numerical differentiation 
is not a significant factor in our calculation. Hereafter, 
we adopt F3, which is expected to cause the least error 
among the three methods. 

Next, we compare the luminosity evolution by varying 
the particle number TV in Fig. [4] As can be seen, how- 
ever, the two peaks are almost independent of the input 
particle number, while the steady emission after and be- 
fore collision shows substantial increases for smaller N. 
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FIG. 4: Dependence on particle number (5.4 x 10 5 : solid, 
5.4 x 10 4 : dashed, 5.4 x 10 3 : dotted). Other parameters are 
the same as the model-0. 
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FIG. 6: Time evolution of luminosity of gravitational waves 
for model-0 (solid), model-A (dashed), and model-D (datted). 




FIG. 5: Time evolution of luminosity of gravitational waves 
for model-0 with contribution from halo (dashed), disk (dot- 
ted), bulge (dash-dotted) only. Total luminosity is also shown 
(solid). 



Thus, we can expect that the feature of the peaks and 
total emitted energy can be studied reliably with particle 
numbers larger than ~ 10 4 . This particle number should 
be compared with that of the simulations in the next sub- 
section, N = 2.4 x 10 5 . On the other hand, the steady 
emission is affected by the resolution of simulation and 
would be highly suppressed in real galaxy merger. 

Having established the above two elements of this 
model, we next study the relative importance of galaxy 
components. Although the mass of a galaxy is mostly 
due to the contribution from its halo, disk and bulge may 
also considerably contribute to gravitational wave emis- 
sion because they are more concentrated than the halo 
itself. To account for this, we calculate the quadrupole 
Jjj, Eq. for each component. However, the dynam- 
ics of particles in all components was solved taking the 
gravitational potential of all the particles into account. 
Fig. [5] shows the contribution of each component (halo, 
disk and bulge). We can see that, at the peaks, the lu- 
minosity is dominated by the contribution from the halo 



and other contributions arc at most 12%. Remembering 
that the mass of the disk is about one seventh of that of 
the halo itself, and that the luminosity is apparently pro- 
portional to the squared mass. The contribution of the 
disk (~ 10%) is larger than a naive expectation. This 
would be due to the fact that the disk is more concen- 
trated than the halo, so that the disk emits gravitational 
waves more effectively. We tested the effect of concentra- 
tion of the disk component by changing the mass profile 
artificially. As a result, a disk with size 80% and 120% 
of the fiducial size contributes 14% and 9% of the total 
luminosity, respectively, while it is 12% for the fiducial 
mass profile. Thus, it is expected that the gravitational 
wave luminosity is dominated by the halo component for 
a reasonable range of the disk concentration. 

Furthermore, we note that the ratio of the luminous 
mass to the dark mass Md a rk/-^ium that we adopted is 
about 7, which is smaller than that of typical galaxies 
(~ 10). Thus, our calculations are overestimating the 
disk and bulge contributions. This fact justifies our later 
calculations with galaxies obtained from a cosmological 
simulation which have only halo components. 

Finally, we compare the gravitational-wave luminosi- 
ties for mergers of modcl-0, model- A and modcl-D 
galaxies. The peak luminosity for modcl-A is 4.7 x 

10 29 erg/sec, which is about one third of that of model- 
0, even though both models have almost the same mass. 
This is because model-0 has a more concentrated halo 
than model-A as can be seen from Fig. [7] On the 
other hand, the peak luminosity for model-D is 2.4 x 

10 30 erg/sec, which is approximately 5 times larger than 
that of model-A. Since the bulge and disk components 
of model-A and D are almost identical, the difference of 
peak luminosities comes from the halo mass components 
of these models. The halo mass distributions of the two 
models are similar within 10 kpc, while model-D has more 
extended halo than model-A as is shown in the bottom 
panel of Fig. 10. Accordingly, the ratio of halo masses 
of model-D and A is about 7.6. Thus, gravitational-wave 
luminosity is very sensitive to both the mass and the den- 
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FIG. 7: Top: the density profile of model-0 (solid), model-A 
(dashed), and model-D (dotted). Bottom: the mass of model- 
(solid), model-A (dashed), and model-D (dotted) within the 
radius. 



TABLE II: Models used in this subsection. Column(l):model 
name. Column(2):mass of one galaxy. Cohmm(3):initial rel- 
ative velocity (relative velocity in the z-direction is set to 
zero). Column(4):initial separation between the galaxies in 
the ^-direction (initial separations in the y- and z-directions 
are set to zero). 

model M(/t~ 1 M Q ) (AV x ,AV y ) (kms' 1 ) Ax (h~ L kpc) 



10 



A-l 


3.1 


i x 10 12 


(0,0) 


1200 


A-2 


3.1 


S x 10 12 


(220,0) 


1200 


A-3 


3.1 


i x 10 12 


(380,0) 


1800 


B-l 


3.1 


S x 10 12 


(0,0) 


850 


B-2 


31 


i x 10 12 


(0,70) 


850 


B-3 


3.1 


S x 10 12 


(0,140) 


850 
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i x 10 11 
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Top: Density profile of the fiducial halo with 1.3 x 
lO^/i^M© (solid) which is well fitted by the NFW profile 
with r B ~ 170 h~ kpc (dashed). For reference, the density 
profile of the galaxy (model-0) used in section IIII Al is also 
shown (dotted). The density and the radius are normalized 



with the average density /9 a v g , and the virial radius 



of 



each galaxy. However, the virial radius of the galaxy used 
in section IIII Al is replaced by the tidal radius. Middle: the 
mass of the fiducial galaxy (solid) and one used in section llll Al 
(dotted) within the radius. As in the top figure the radius is 
normalized, and the axis of ordinate is normalized by the total 
mass of each model. Bottom: rotation curves of the fiducial 
galaxy (solid) and one used in section [III Al (dotted) . 



B. Simulations with Type B Initial Conditions 



sity structure of halo, making quite clear the importance 
of using a realistic model. In the next subsection, we sim- 
ulate mergers of galaxies with a NFW-like density profile 
and study the mass dependence. 



In this subsection, we adopt a dark halo from a cos- 
mological simulation (for details, see [HI) to prepare ini- 
tial conditions. The halo identification is done in a two- 
step manner. First, we select candidate objects using 



6 



the friends-of- friends (FOF) algorithm [l4j]. We set the 
linking parameter b = 0.2. Secondly, we apply the spher- 
ical over density algorithm to the located FOF groups. 
To each FOF group we assign a mass such that the en- 
closed mass within the virial radius is A x p CI it(z), where 
Pcrit(z) is the critical density for which the spatial geom- 
etry is flat. Based on the spherical collapse model fl5j |. 
A is 200. The adopted halo has mass 1.3 x 10 14 /i -1 M Q , 
particle number ~ 2.4 x 10 5 , and a density profile well 
approximated by the NFW profile, 



P(r) 



{r/r s )(l+r/r s y 



(4) 



where r s ~ 170 /i -1 kpc (Fig. [5]). We placed two identi- 
cal halos with various initial relative velocities. Also we 
vary mass of one or both of two galaxies to study the 
dependence of mass (ratio) by a scaling law based on an 
empirical relation obtained in [lj|, M V - 1T oc V„" ax , where 
a ~ 3.4. Here L, V max and M v ; r are luminosity, maxi- 
mum circular velocity and virial mass, respectively. This 
relation is realistic taking into account galaxy rotation. 
From this relation, the following relation is obtained [l7| , 



M 
Mo 



>\ 1-2/ a 



(5) 



where ro and Mq are the original radius and mass, and 
r' and M 1 arc the scaled radius and mass, respectively. 

Table [TT] summarizes the models used for our calcula- 
tions. Models A-l, A-2 and A-3 are used to study the 
dependence of initial relative velocity in the direction of 
initial separation (x-direction). Effects of initial relative 
velocity in the y-direction, that is, the relative angular 
momentum, are studied by models B-l, B-2 and B-3. Fi- 
nally, models C-l and C-2 have different mass ratio, 1/10 
and 1/100, respectively. 

Fig. [3] shows the comparison of luminosity evolution 
varying the initial relative velocity in the direction of ini- 
tial separation (models A-l, A-2 and A-3). Here, the 
maximum relative velocity with which two galaxies are 
gravitationally bound is about 220 km sec -1 . There- 
fore, model A-2 represents a case where two galaxies 
are marginally bound, while model A-3 represents a case 
where two galaxies are not bound and pass through each 
other after the first collision. This is why there is only 
one peak in the luminosity of model A-3. Other than this 
point, we cannot see any differences between the three 
models in this figure that is, the peak luminosity and 
width look almost independent of AV X . 

To see the difference between these three models more 
closely, we plot the amplitude of the plus mode in the 
z-direction in Fig. [TD] The amplitude is calculated at 
10 /i -1 Mpc from the galaxies and the peak amplitude is 
of order 10 -13 . The difference in the amplitudes is about 
20 — 50% and increases as AV X increases. 

Fig. [H shows the amplitude of the plus and cross 
modes in the x-, y- and z-dircctions for model A-2. We 
see that gravitational waves are mostly emitted as the 




FIG. 9: Comparison of luminosity evolution varying the initial 
relative velocity in the direction of initial separation (models 
A-l, A-2 and A-3). 
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FIG. 10: Amplitude of the plus mode in the z-direction at 
10 /i -1 Mpc from the galaxies for models A-l, A-2 and A-3. 



plus mode in y- and z-directions, and the cross mode is 
smaller by one order. This comes from the fact that the 
galaxy motion is almost axisymmetric with respect to the 
x-axis. The small difference between radiation of the y- 
and z-directions is due to the small deviation of galaxy 
structure from the spherical symmetry. 

In Fig. 1121 we show the comparison of luminosity evo- 
lution varying the initial relative velocity in the direc- 
tion orthogonal to the direction of the initial separation 
(models B-l, B-2 and B-3). When the two galaxies have 
relative angular momentum, it takes more time to col- 
lide and relax into a single galaxy. Thus, the peak width 
broadens for large AV y while the peak luminosity is sup- 
pressed. 

Fig. Q2] is the spectrum of the luminosity for models 
B-l, B-2 and B-3. Here the spectrum is calculated by the 
discrete Fourier transform, 



JV-l 



£gw(/«) = At 2J LGw(*fe)exp 



k=0 



2-Kikn 
N 



(6) 



where n = 0, • ■ ■ , N — 1 and N is the number of time 
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FIG. 11: Amplitude of the plus (top) and cross (bottom) 
modes in the x-, y- and z-directions at 10 /i _1 Mpc from the 
galaxies for model A-2. 




FIG. 12: Luminosity evolution varying the initial relative ve- 
locity in the direction orthogonal to the direction of the initial 
separation (models B-l, B-2 and B-3). 
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FIG. 13: Fourier spectrum of the luminosity for Models B-l, 
B-2, B-3. 
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FIG. 14: Amplitude of the plus (top) and cross (bottom) 
modes in the x-, y- and z-directions at 10 h~ Mpc from the 
galaxies for model B-2. 



bins. Here At = 5/h Myr is the width of each time bin 
and tk — kAt is the time of the fc-th bin. Gravitational 
wave frequency is represented by /„ = n/(NAt). We 
can see that increasing AV y decreases the characteristic 
frequency which corresponds to the peak width, as is 
expected from Fig. [12] 

Relative angular momentum also affects the direction 
dependence of the emission. Fig. H4l shows the amplitude 
of the plus and cross modes in the x-, y- and z-dircctions 
for model B-2, which correspond to Fig. [TT]for model A- 
2. We see that the plus and cross modes are emitted with 
the same order of magnitude in the z-direction, while 
there is no significant difference between model A-2 and 
B-2 in the x- and y-dircction. This is because the galaxy 
motion is no longer axisymmetric, due to the relative 
angular momentum. 

Fig. [T5] shows the dependence of the peak luminosity 
as a function of the galaxy mass M fixing the mass ratio 
to 1 : 1. As can be seen, the luminosity is well fitted by 
M 10 / 3 - 4 . One may expect that the luminosity is propor- 
tional to M 2 because the quadrupolc is proportional to 



M. However, this is not the case because the length and 
time scales are also dependent on M. We can understand 
this behavior as follows. 

Defining the typical length and time scales as L and T, 
respectively, the quadrupolc and luminosity are roughly 
given by, 



ML 



ML 2 



Lgw 



M 2 L 4 

J^6 



(7) 



On the other hand, assuming the gravitational potential 
energy between the two galaxies is equal to the kinetic 
energy of their relative motion at the collision, charac- 
teristic relative velocity V ~ L/T is evaluated as, 



M 

T' 



Thus, we have, 



L 3 
M' J 



(8) 



(9) 
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FIG. 15: Peak luminosity as a function of mass fixing mass 
ratio to 1 : 1. 
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"GalactlCS" which consist of disk, bulge and halo. Wc 
showed that the features of peaks such as the width and 
luminosity are reliably simulated with particle numbers 
larger than ~ 10 4 . Following this, the relative impor- 
tance of components was investigated and it was found 
that the dominant contribution to gravitational- wave lu- 
minosity comes from halo component, while contribu- 
tions from other components were less than 10%. 



Second, we used initial conditions with a dark halo 
adopted from a cosmological simulation. This halo had 
a realistic density structure and was more suitable for a 
precise estimation of gravitational-wave luminosity. The 
peak luminosity amounted to 10 31 erg sec -1 for the col- 
lision of two halos with mass 3.8 x 10 12 ft, _1 M . Wc 
showed that the initial relative velocity in the direction of 
the initial separation does not significantly affect gravita- 
tional wave emission, with a difference of about 20 — 50% 
in amplitude. On the other hand, we found that the ini- 
tial relative angular momentum broadens the peak width 
and suppresses the luminosity by one order. Mass de- 
pendence of the peak luminosity was also investigated. 
In contrast to the naive expectation, Lew °t M 2 , we 
obtained Lgw oc M 3 . We gave a simple analytic inter- 
pretation of this behavior based on the scaling relation 
of the mass and size of galaxies. Also this analysis was 
shown to be applicable to the dependence on the mass 
ratio. 



FIG. 16: Luminosity evolution varying the mass ratio (models 
B-l, C-l and C-2). 



and then, 



From Eq. 



M 2 L 4 

y6 



M 5 



(10) 



), we have L oc M 1 2 l a and finally obtain, 
L Gw cxM^, (11) 



with a ~ 3.4. 

Fig. [TB] shows the comparison of luminosity evolu- 
tion varying the mass ratio. The peak luminosity is 
roughly proportional to (M1/M2) 1 " 5 where Mi < M2. 
This dependence may be understood by replacing M with 
\JM\M2 in the discussion in the previous paragraph. 



IV. SUMMARY AND DISCUSSION 

In this paper, we performed a systematic study of grav- 
itational waves from galaxy mergers through TV-body 
simulations with Gadget2. Two types of initial condi- 
tion were adopted. First we used galaxies generated by 



To argue the observability of gravitational waves from 
galaxy merger through the B-modc polarization of CMB, 
it is necessary to evaluate gravitational wave background, 
which is contributed to by the entire history of galaxy 
merger. In future research we will present a detailed es- 
timation by combining the elementary process studied 
here and the mer ger history based on hierarchical galaxy 
formation model [l8j. Here we give a very rough order- 
of-magnitude estimation based on several simplifying as- 
sumptions. Let us assume that all the current galaxies 
have mass of 10 12 /i -1 M© and were created through suc- 
cessive mergers starting from the progenitor galaxies with 
mass 10 lo /i _1 A/ Q . Further, we assume that all the merg- 
ers took place with mass ratio 1:1 and head-on orbit, 
and neglect the redshift of gravitational waves. Count- 
ing the number of mergers in the current horizon volume 
and multiplying the gravitational-wave energy emitted 
at each merger taking the scaling law in Fig. [IS] into ac- 
count, we obtain the density parameter of gravitational 
wave background as f^cw = Pgw/Pc ~ 10~ 21 where 
p c is the critical density. This corresponds to inflation- 
ary gravitational waves with the tensor-to-scalar ratio 
~ 10~ 6 . Although it would be too small to be detected 
through B-mode polarization of CMB by the existing in- 
struments, delensing of CMB polarization may allow us 
to probe such weak gravitational wave background in the 

future mm. 



9 



Acknowledgments 

TI appreciates Vicent Quilis, and A. Cesar Gonzalez- 
Garcfa for the detailed information on their simulation 
Q. This work is supported in part by JSPS Grant- 
in-Aid for the Global COE programs, "Quest for Fun- 
damental Principles in the Universe: from Particles to 



the Solar System and the Cosmos" at Nagoya University. 
KT is supported by Grant-in- Aid for Scientific Research 
No. 21840028. NS is supported by Grant-in-Aid for Sci- 
entific Research No. 22340056 and 18072004. This re- 
search has also been supported in part by World Premier 
International Research Center Initiative, MEXT, Japan. 



[1] S. Mollerach, D. Harari, and S. Matarrese, Phys. Rev. D 
69, 063002 (2004), arXiv:astro-ph/0310711. 

[2] K. N. Ananda, C. Clarkson, and D. Wands, Phys. Rev. 
D 75, 123518 (2007), arXiv:gr-qc/0612013. 

[3] D. Baumann, P. Steinhardt, K. Takahashi, and K. Ichiki, 
Phys. Rev. D 76, 084019 (2007), arXiv:hep-th/0703290. 

[4] C. Carbone, C. Baccigalupi, and S. Matarrese, Phys. 
Rev. D 73, 063503 (2006), arXiv:astro-ph/0509680. 

[5] C. Lacey and S. Cole, MNRAS 262, 627 (1993). 

[6] C. Lacey and S. Cole, MNRAS 271, 676 (1994), 
arXiv:astro-ph/9402069. 

[7] V. Quilis, A. C. Gonzalez-Garcfa, D. Saez, and J. A. 
Font, Phys. Rev. D 75, 104008 (2007), 0704.3009. 

[8] V. Springel, N. Yoshida, and S. D. M. White, New As- 
tronomy 6, 79 (2001), arXiv:astro-ph/0003162. 

[9] V. Springel, MNRAS 364, 1105 (2005), arXiv:astro- 
ph/0505010. 

[10] K. Kuijken and J. Dubinski, MNRAS 277, 1341 (1995). 
[11] W. Jaffe, MNRAS 202, 995 (1983). 

[12] J. F. Navarro, C. S. Frenk, and S. D. M. White, Astro- 



phys. J. 462, 563 (1996), arXiv:astro-ph/9508025. 
[13] S. Masaki, N. Yoshida, and M. Fukugita, in preparation 
(2010). 

[14] M. Davis, G. Efstathiou, C. S. Frenk, and S. D. M. White, 

Astrophys. J. 292, 371 (1985). 
[15] P. J. E. Peebles, The large-scale structure of the universe 

(1980). 

[16] J. S. Bullock, T. S. Kolatt, Y. Sigad, R. S. Somerville, 
A. V. Kravtsov, A. A. Klypin, J. R. Primack, and 
A. Dekel, MNRAS 321, 559 (2001), arXiv:astro- 
ph/9908159. 

[17] A. C. Gonzalez-Garcfa and M. Balcells, MNRAS 357, 

753 (2005), arXiv:astro-ph/0411645. 
[18] M. Enoki, K. T. Inoue, M. Nagashima, and N. Sugiyama, 

Astrophys. J. 615, 19 (2004), arXiv:astro-ph/0404389. 
[19] U. Seljak and C. M. Hirata, Phys. Rev. D 69, 043005 

(2004), arXiv:astro-ph/0310163. 
[20] K. Sigurdson and A. Cooray, Physical Review Letters 

95, 211303 (2005), arXiv:astro-ph/0502549. 



